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Abstract 

We present the one-loop corrected decay widths for the decays of the neutral 
Higgs bosons h°, H° and ^4° into a neutralino pair %%i Xn ( m > n = 1, . . . , 4) and to the 
decay Xm ~~ *■ (^°> H > ) + Xn- The corrections contain the onedoop contributions 
of all fermions and sfermions. All parameters are taken on-shell. This requires a 
proper treatment of the neutralino mass and mixing matrix. The dependence on 
the SUSY parameters is discussed. The corrections can be large in certain regions 
of the parameter space. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) [|IJ is considered the most attrac- 
tive extension of the Standard Model. The MSSM requires the existence of two isodoublets 
of scalar Higgs fields, implying three neutral Higgs bosons, two CP-even bosons (h°, H°), 
one CP-odd (A ), and two charged Higgs bosons {H^. Searching for these Higgs bosons 
is one of the main goals of all future colliders as the Tevatron, LHC, and an e + e~ linear 
collider. The search strategies very much depend on the way these Higgs bosons decay. 
It is therefore mandatory to have a clear picture of the decay modes. Thus it is necessary 
to calculate the widths and branching ratios of the various decays as precisely as possible. 

The lightest Higgs boson h° with a mass of at most 140 GeV will decay mainly into bb 
and to a lesser extent into t + t~. It is, however, possible that it also decays as 

h ^x°, + xl (1) 

where Xi is the lightest neutralino. In the case of .R-parity conservation, this decay is 
invisible, and its appearance would reduce the branching ratios of the other decay modes. 
The heavier neutral Higgs bosons H° and ^4° may decay into a pair of neutralinos 

(H°,A )^x m + X° n , (2) 

with (m,n = 1, ... ,4). At tree level, the decays occur by higgsino-gaugino-Higgs boson 
couplings 0, and are therefore sensitive to the components of neutralinos. The decays 
(H) and (H) as well as those of H ± — > xtXmi = 1>2) have been numerically analyzed 
in H El at tree level. Electroweak corrections to the widths of H ± — » xtXm due to 
one-loop exchanges of the third generation quarks and squarks were recently calculated 
in H. The one-loop corrections, involving fermions and sfermions, to the invisible width 
of (h°, H°, A ) — > Xi + X? have been calculated in the higgsino limit of xi, (H < Mi,M 2 ) 
in ||, and in the gaugino limit of Xi, (I A 4 1 ^> Ml, Ma) very recently in f7|. (Here Mi and 
M 2 are the £7(1) and SU(2) gaugino mass parameters, respectively, and \i is the higgsino 
mass parameter.) In these limiting cases, the wave-function corrections can be neglected 
and no renormalization is necessary. The couplings of (h°, H°, A ) to xi Xi a ^ so enter 
in the neutralino-quark interaction 0, a process which is very important for the dark 
matter search || [J, where one looks for the elastic scattering of neutralinos x,i °ff nuclei 
in a detector. Moreover, since the decays (|l],|D are generated by gaugino-higgsino-Higgs 
boson couplings at tree level, they can be also useful to probe the components of the 
neutralinos, complementary to the pair production process e + e~— > XmXn EH- 

In this paper, we present the one-loop corrections to the widths of the decays (P and @ 
due to the exchange of all fermions (quarks and leptons) and their superpartners (sfermions) 

The decays (P and @ are particularly interesting because the calculation of their ra- 
diative corrections requires corrections to the neutralino mass matrix and mixing matrix 
in addition to the conventional wave-function and vertex corrections with counter terms. 
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The one-loop corrections to the neutralino mass and mixing matrix in the on-shell renor- 



malization scheme were already worked out in |Tl[ and they will be used here. 
Related to these decays are the decays of neutralinos into Higgs bosons, 



X° 

Am 



(h°,H°,A° 



+ X° n 



(3) 



These decays are also important as they occur in the cascade decays of gluinos and/or 
squarks, g — > qqx° m an d Q with Xm then decaying according to (13). The decays 



(0) with a real Higgs boson emission |TJ, [L3| as well as three-body decays due to an off- 
shell Higgs boson |[3| have been studied at tree level. In this paper, we also present the 
formulae for the decays (^j including the one-loop corrections. 



2 Tree-level widths 

Throughout this paper, we will use the notations m^o = rrii and H k = {h°, H°, A , G }. 
In a non-unitary gauge we have the ghost G°. The momenta are assigned as (k = 1, 2, 3; 
m, n = 1, . . . , 4) 

H° k ( P ) ^ X°m(h) + X°n(k2) ■ (4) 

All couplings are given in the Appendix [A] (or it is referred to previous works). 
The tree-level widths for a neutral Higgs decaying into two neutralinos is |3J 

■ptree ptree / itO . ~ o ~ o \ 

a 2 

K{rn H o,m m ,m n ) \F mnk \ [m H o -m m -m n - 2(-l) ks m m m n \ , (5) 



87im 3 H0 (1 + 5 mn ) 



k 



with k(x, y, z) = ((x — y — z) 2 — Ayz) 1 / 2 . The couplings F^ nk are given in the Appendix, 
eqs. (OS)- 

For the decay of a neutralino into a lighter one and H k , we get 



-ptrec _ -ptree/ ~0 , ttO ~0\ _ 
1 X° ~ L ~^ n k Xn) — 

16?rm 3 Wmnk? [^m ~ m \ + m l + 2(-l) 5k3 m m m n \ . (6) 

In our convention, the 4x4 neutralino mixing matrix Z, which diagonalizes the neutralino 
mass matrix Y, is real. Therefore, the neutralino mass parameters m m and m n can be 
positive or negative. 
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3 One-loop corrections 



We calculate the one-loop corrections to the amplitudes of the decays @ stemming from 
fermion and sfermion exchange. The renormalization is done in the on-shell scheme. All 
one-, two-, and three-point functions [|15| used for calculating the loop integrals are given 

m. 



in the convention 



The correction to the coupling F^ nk is 

mnk mnk ' mnk i \ ' ) 

with the ultraviolet (UV) finite one-loop correction 

^ F nmk = N c ( SF mnk + SF mnk + SF mnk) > ( 8 ) 

flavors 

with the color factor N[ = 1 for (s)lepton and N[ = 3 for (s)quark exchange, ^flavors 
stands for the summation over all (s)fermion flavors, e. g. (top, stops), (bottom, sbottoms), 
(tau, staus), etc.. For convenience, the color factor N[ is given only in the total correction 
term eq. 

In our convention, both F^ nk and AF^ nk are real. Therefore, the corrected widths can 
be written as 

/ A PO \ 2 

■pcorr. -ptrcc / i , mnk \ /r\\ 

L p - L p l i + _ ro I ' W 

\ mnk / 

with the decaying particle p = H k or 



The vertex correction stems from the two diagrams shown in Figs. |T]a and [3]b. Because 
of the Majorana nature of the neutralinos the charge conjugated (s)fermion fields denoted 
by the superscript "c" can also circulate in the loop. 
For h° and H° (a = 1, 2) we have 



si 2 



+ 



9 SF mna = J^y H | m / (4m4n + Hm h in) ( m m + Wn)^ + l{m m C\ + m n C, 

+ ( a imH n + a{rMm) ( m / + m \ + m m m n ^C l + (m m + m n ) (m m C\ + m n C l ^j + B 



1 ±ci 



(4tt) 2 



i,j =1 

For A the vertex correction reads 
isl 2 ' 



(a{ m bj n + ajJji^rrifCj — {a{ m a^ n + b{ m b^ (m m C^ + m n C^ 



(10) 



9 $ F !nL = 7^T2 £ | m / (pLflL + HmHn) ( m m + ™n) ^0 + 
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+ ( a im b in + a LHm) ("^m"^n ~ ™| + TJlfj Cq + [ffl„ - 77i m ) (m m CJ - m n C. 

+77 \ 2 ( a i m bjn ~ a jntfm) m f ^0 + ( a im a jn ~ HmPjn) {^mPl ~ m n C^ 

l 47r J ij=l L 



fl (m*o,mj,mj), 



The abbreviations Bq 
The wave-function correction is given by 



C 



C.. (m^, m^o , m 2 ,, m| ,m 2 f ,m 2 f ) 



III) 



and 



have been used. 



0(io) 



mnk 



A 7^° P° -L A 7* F° -L A 7* TT" 

Ik Tfinl * nm nn.k ' ^ an m 



qm qnk 



J qn mqk 



(12) 



with the implicit summations over I = 1,2 for k = 1 or 2, / = 3,4 for k = 3, and g = 
1, . . . , 4. 5Z* are the wave-function constant terms for the neutralinos given in (|20|), (pT|). 
The wave-function constant terms for the Higgs bosons 5Z H ° are 



SZ, 



SZ, 



kk 

H° 
kl 



2 2 

TflrrO TTlrrO 



^\K\m 2 H? )-K\m 2 H0 ) 



l^k, 



(13) 
(14) 



with k,l = 1,2 for the system #°) and k,l = 3,4 for Eq. (ID has been 

symmetrized with respect to (k, I). This is due to the on-shell renormalization of the 
Higgs mixing angle a (k,l = 1,2) or (3 (k,l = 3,4). In this scheme (IP-HI, extending |17] 
for quark and lepton mixing) the counter terms for the mixing angles are determined by 
the requirement that they cancel the antisymmetric parts of the wave-function corrections. 

The decays of A are a little complicated by the contribution of the A — ZJ° mixing 
in addition to the A — G° mixing in eq. (0). Moreover, both depend on the gauge 
parameter £. However, the sum of these two contributions is independent of £, as it is 
shown in Appendix C. Here we work in the £ = (Landau) gauge, where the contribution 
of the A — Z° mixing vanishes, and use (|T4|) with m H o = 0. The resulting on-shell tan/? 
agrees with the one defined by the A°-Z° mixing [ 18j [19], |Il 



The Higgs self-energy contributions due to fermions and sfermions are written as 



+ n 



kl 



(k 



kl 



The fermion contribution Yl H ° ^ (Fig. [TJc) is 



K b w (k 2 ) 



n 



if (a) 



cci 



(^ 



(4tt)' 

(47T) 2 



(fc 2 - 4mj) B (k 2 , m}, m)) - 2A (m 2 f ) 



k 2 B (k 2 , m), m)) - 2A (m 2 f ] 



(15) 

(16) 
(17) 
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Figure 1: One-loop Feynman graphs with fermion and sfermion exchange contributing 
to the neutral Higgs boson-neutralino-neutralino decay amplitude. The superscript "c" 
denotes the charge conjugated states. 
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with a,b — 1,2 and c,d = 3,4. The sfermion contributions U H ° ^ (Fig. |l|d) and U H ° ^ 
(Fig. |l|e) are 



n 



s 0( v 



(4vr) 2 iJ=1 



2 ^ 



n 



H° (c) 
kl 



(47T) 



i=l 



(19) 



where fc, / = 1,2 or 3,4. Tjy in eq. ( |ToD represent momentum-independent contributions 
from the tadpole shifts [0, pl| and leading higher-order corrections. We include the latter 
by the renormalization group improvement as in Ref. ||20l. Since the zero-momentum 



contribution n^°(0), including T^i, is very large it is often resummed as in Refs. |H| |2T |. 
In practice, we calculate the effective m#o(a = 1,2) and a obtained from the effective mass 
matrix, which includes the IF^f \k 2 = 0) contribution with m^o, tan/?, and the (s)quark 
parameters, and regard them as the lowest-order parameters. If one is replacing a in all 
the previous formulae with the effective one, the self energies IT^ (k 2 ) in the wave-function 
correction and 5a must be replaced by AI1^ (A; 2 ) = IT^ (k 2 ) — n^°(0). Nevertheless, the 
form of their sums eqs. ( |I1^|T4"1) is not affected by the elimination of 11^ (0). 



The neutralino wave-function terms read 



^1 



- Re <j IT* p (m p ) + 2m p m p Tl* p {m p ) + m/EC 
1 



m p — m s 



Re< m 



n ps(^J - n p S (^ P ) + m f 



n|/(m; 



) - n* s 

I ps 



ps \ m pJ 



(20) 
(21) 



p s. As before, 5Z* S in ( pT|) has been symmetrized by subtracting the counter term 
for the rotation matrix Z of the neutralinos [|ll|] . The neutralino self-energies due to the 
sfermion-fermion loop (Fig. |l]f ) are 

E(4°£ + 6 W) m l m \) > ( 22 ) 



n? s (A; 2 ) 



1 



(4 



7T 



t=l 



n*/(fc 2 ) 



(4vr) 2 



E(4A4 + 4s h ip) Bo(k 2 , mj, m|) 



(23) 



i=l 



see also pCj. 



We need the counter term for the couplings gF^ nk , which is a function of the gauge 
couplings g, g', the Higgs boson mixing angle a (for k — 1,2) or /3 (for k — 3), and the 
neutralino rotation matrix Z, as shown in eq. ( |A.2| ). The counter terms for (a, (3) and Z 
in the on-shell scheme [[ll|] are already included in the wave-function corrections ([14]) and 
respectively. The remaining counter term of 5g and 5g' is, after being absorbed into 



the correction to F^ nk , 



SR 



0(c) 



rank 



E d% 



x=3,4 



mx Zj n2 



7 V \ 8 A 
9 



y Z mx Z n i + Z m iZ nx ) tan 8w — r 

9 



.(24) 
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We fix the electroweak gauge boson sector by mz, mw, and e. One gets from the relations 



g = e/%, g' = e/c w , and c w = m w /m z (c w = cos6 w , s w = sm6 w ) |22], [16 ] 



Sg _6e c 2 w { 5m 2 w 5m z \ 5g' _ 5e 1 / 5m 2 w 5m z \ 

g e 2 sly \ m 2 w m| / ' g' e 2\m 2 v m\ ) ' 

The formulae for 5m\y and Smz can be also found in [|TT|] and for Se/e in the Appendix |B|. 

Now all parts are given which are needed in order to calculate the (UV finite) one- 
loop contribution to the neutral Higgs boson-neutralino-neutralino coupling, eq. (|8|). The 
vertex correction part SF^l is given by eqs. (|T0"1) and (0), the wave-function correction 
term 6F%$ by eq. (pD, and SF^. by eq. (||). 

Further, one has to note that the on-shell masses and the mixing of the neutralinos are 
not independent of each other. In fact, when the gauge and Higgs boson sectors are 
fixed, the neutralino sector is determined by three free parameters only. Here we follow 
the method given by [[II]]: The on-shell mass parameters M and fi are defined as the 



elements of the on-shell mass matrix X of charginos, and the on-shell mass parameter 
M' is defined as the element of the on-shell mass matrix Y of neutralinos . The finite 
correction AY = Y — Y trcc , where y tree is the tree- level mass matrix in terms of the 



on-shell parameters (M, /x, M', mz, sin^, tan/5), is calculated by eqs. (42-51) in JTT 
The one-loop corrected on-shell masses rrii and mixing matrix Z = Z tTee + AZ are then 
obtained by diagonalizing Y. For a proper treatment of the loop corrections, the resulting 
shifts of the masses and the mixing matrix from the tree-level values have to be taken 
into account. 



4 Numerical results 

For simplicity, we will take in the following (if not specified otherwise) for the soft breaking 
sfermion mass parameters of the first, second and third generation Mq i = Mfj 2 = 

M D,a = = = M Qs = f M u 3 = Ti M n 3 = M U = M e 3 = M Q = 500 GeV 

and for the trilinear couplings A t = Ab = A T = A = 300 GeV. We take mz = 91.2 GeV, 

m w = 80 GeV, m^o = 500 GeV, m t = 175 GeV, m h = 5 GeV, and m T = 1.8 GeV. Masses 

of all other SM fermions are neglected. We use the GUT relations M' = | tan 2 8w M 

and for the gluino mass rrig = {a,s/ a Ew) sin 2 @w M. The other input parameters are 

{tan/3, M, //} (all as on-shell parameters). For the values of the Yukawa couplings of 

the 3 rd quark generation (h t , hb), we take the running ones at the scale of the decaying 

particle mass. 

In our numerical analysis we have discussed four cases: the tree-level width, the corrections 
(0 H) with the tree-level Z and m 8 ("conventional correction"), the corrections (|7|~fl) with 
the one-loop corrected Z and tree- level m ; ("conventional + AZ correction"), and the 
corrections (0-^) with the one-loop corrected Z and one-loop corrected rrii (full correction). 



S 



The "conventional correction" corresponds to the correction to the gaugino-higgsino-Higgs 
boson coupling, "conventional + AZ correction" includes the correction to the neutralino 
components, and the correction due to the shift of rrii is added in the full correction. 



In Fig. 0a we show, as a function of fi, the tree-level widths of H° — > X1 + X25 H° X1 + X3 



Xi + xi 

and H° — > X2 + X%i respectively, for tan/3 = 10 and M = 150 GeV. The H° mass is 
rrifjo ~ 500 GeV. The widths vary with the gaugino and higgsino components of the 
various neutralino states. Fig. 0b exhibits the corrections to the width of H° — > Xi + 
The "conventional", "conventional +AZ", and full corrections are shown. One can see 
that, compared to the "conventional" correction, the corrections by the shifts AZ and rrii 
cannot be neglected. Figs. 0c and §d show the corrections to the widths of H° — > Xi + X3 
and if — > x\ + x!L respectively. While the "conventional" correction is dominant for 
\i < in Fig. 0c, the "AZ" correction is dominant in Fig. |2|d. The full corrections amount 
to several %. 
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Figure 2: Tree-level widths (a) of the decays H° — > V1 + X2 (solid), H° X1 + X3 (dashed) 
and H° — > X2 + X2 (dotted) and corrections to the widths of these decays (b), (c), and (d), 
respectively, as a function of /1 for tan ft = 10 and M = 150 GeV. The full, dashed, dash- 
dotted line corresponds to the full, "conventional + AZ" , and "conventional" correction. 
The grey areas are excluded by the bounds m~± > 100 GeV, m h o > 92 GeV. 
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In Fig. [| we show the tree and corrected widths of (a) H° — > X1+X2 decay with tan (3 = 10, 
M = 500 GeV and /i = 150 GeV, and those of (b) A -> xl + X° 3 decay with tan (3 = 50 
and M = /1 = 300 GeV, as functions of Mq. In Fig. |a, the decay is suppressed due 
to the small gaugino components of Xi an d xH- The "conventional+AZ" correction is 
close to the full correction and therefore not shown here. We see that the "conventional" 
correction is dominant. In contrast, the "Am" correction in Fig. ^0 is large and negative 
(up to —16%), which dominates over the positive "conventional" correction (up to +4%). 
This is because the decay in Fig. |3|b is kinematically suppressed and sensitive to the shift 
of m r We note that the sfermion loop corrections do not decouple in large Mq limit, 
due to the supersymmetry breaking corrections [2~3~| to the gaugino-higgsino-Higgs boson 
couplings. 




200 300 400 500 

Mq [GeV] 



600 300 400 500 

Mq [GeV] 



600 



Figure 3: The widths of the decays H° — > xi + X2 ( a ) an d ^4° —* Xi + x\ (b) as a function 
of Mq. The dotted line corresponds to the tree-level width, the dash-dotted, dashed, and 
solid line corresponds to the "conventional", "conventional + AZ n , and full correction, 
respectively. The parameters are tan/5 — 10, M — 500 GeV, and \i = 150 GeV (a) and 
tan/3 = 50 and M = /i = 300 GeV (b). 

Figs. |]a and f|b show the dependence of the widths on the trilinear coupling A for the same 
decays modes and parameter sets as in Figs. |3]a and |3|b, respectively. The A dependence 
is mainly caused by A t and numerically important in general. 



We also discuss the related decays (Q) of the neutralinos. In Figs. ||a and ||b we show the 
corrections to the width of the decays X3 ~ ¥ h° + Xi an d xl h° + Xi as functions of 
Mq and A, respectively. The parameters are as in Fig. |3|a for Fig. |5|a and tan j3 = 10 and 
M = fi = 300 GeV for Fig. |5|b. The total correction can go up to 25%. 

Finally, Fig. |] shows the \i dependence of the width of the decay h° — > Xi + Xi (HI)) both 
the tree-level value (Fig. |6|a) and the relative one-loop full correction (Fig. |j|b). This 
decay occurs when Xi is sufficiently light and is mainly a U(l) gaugino to escape from the 
present direct search. In order to realize this case, we consider very small M' and take 
the following parameters which are similar to those in Ref. J7|: M Q 



M 



c/1,2 



M 



Di 
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A [GeV] A [GeV] 

Figure 4: The widths of the decays H° — > \\ + X2 ( a ) an d ^° ~> Xi + X3 (b) as a function 
of A. The dotted line corresponds to the tree-level width, the dash-dotted, dashed and 
solid line corresponds to the "conventional", "conventional + AZ" , and full correction, 
respectively. The parameters are tan/3 = 10, M = 500 GeV, and fi = 150 GeV (a) and 
tan/3 = 50 and M = fi = 300 GeV (b). 
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Figure 5: The widths of the decays xt ~ * h° + Xi ( a ) an d Xa ~ y h° + Xi (b) as a function 
of Mq (a) and A (b). The dotted line corresponds to the tree-level width, the dash- 
dotted, dashed and solid line correspond to the "conventional", "conventional + AZ", 
and full corrections, respectively. The parameters are tan (3 — 10, M — 500 GeV, and 
\i = 150 GeV (a) and tan/3 = 10 and M = \i = 300 GeV (b). 

M Q3 = M 0s = M bz = 2M ii2 = 2M£ 12 = 2M U = 2M^ = 500 GeV, A t = 1000 GeV, 
A b = A T = GeV, M' = 30 GeV, M = 120 GeV, m- g = 500 GeV, tan/3 = 20, and 
rriAO = 300 GeV. The loop correction can be comparable to or even larger than the tree- 
level width as observed in Ref. 0. Although the decay width is much smaller than the 
other modes, the effect of the loop correction might be seen in precision studies of h° at 
a linear collider |7| and in neutralino dark matter search [J71 KL B|. 
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/x [GeV] n [GeV] 

Figure 6: The tree- level width (a) and the relative full correction (b) of the decay 
h° — > Xi + Xi as a function of /x for the parameters tan/3 = 20, {M, M', m§, m^o} = 
{120, 30, 500, 300} GeV, {A t , A b , A T } = {1, 0, 0} TeV, and M Qi2 = M Dl 2 = M Di2 = 
M Qa = M 0a = M D3 = 2M Li 2 = 2Mz h2 = 2M ia = 2M^ = 500 GeV. 



5 Conclusions 

We have presented the calculation of the one-loop corrections to the decays 
(h°, H°, A ) ^x° m + x° n and xl ^ {h°, H°, A ) + x° n , (m, n = 1, . . . , 4), with all fermions 
and sfermions in the loop. These decays are special in the sense that they require partic- 
ular care in the treatment of the neutralino mixing and mass matrix in a scheme, where 
all parameters in the neutralino mass matrix Y and mixing matrix Z are defined on-shell. 
We have shown the importance of the corrections to these matrices in addition to the 
conventional corrections (vertex and wave- function corrections with counter terms). We 
have studied the dependence on the parameters M, /i, A, Mq, and tan/3. The corrections 
to the widths of the decays (h°, H°, A ) — > x^L + Xn can S° U P t° ~ 15%, those of the 
decays xl — ► (h°, H°, A ) + Xn to ~ 25%. For the invisible decay h° — > Xi + Xv giving 
up the GUT relation for M', one even gets corrections up to 140%. 
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Appendix 

In the following we give the formulae for the couplings and for 5e. Furthermore, the proof 
of the gauge independence of the processes considered will be given. 



A Coupling parameters 



The if" xl Xn interaction is given by 



C = --qH°Y° F° y° + - o if y° F° 7 5 y° 

2 " a Am mna An 1 <-> i> c Am M mnc I An i 



(A.1) 



with a = 1, 2, c = 3, 4, and 



rank 



+ 



4 



Z m zZ n 2 + Z n3 Z m2 — tan 9y/ (Z m3 Z nl + Z n3 Z ml ) 



(A.2) 



The Z mn are the elements of the neutralino mixing matrix which diagonalizes the neu- 
tralino mass matrix and 



d\ = dk = ( — cos a, —sin a, cos/?, sin/^, 
dk = dt — ( — sin a, cos a, — sin/5, cos P) k - 

The superscript "u" denotes an up-type and "gP a down-type fermion. 

The neutral Higgs boson-fermion-fermion couplings, defined by C = s{ if ° / / are 



(A.3) 



with a = 1, 2 and c = 3, 4, using the Yukawa couplings 

9 m u gm d 



m w sin (3 



h d = 



m w cos P 



-^d d a (A.4) 



(A.5) 
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The fermion-sfermion-neutralino coupling parameters a{ p and b{ p (i = 1, 2, p = 1, . . . , 4) 



have the form 



hp 



gflpRL + hfZ px R f i2 , 

Hp — 9fRpRL + hfZ px R a , 



(A.6) 
(A.7) 



with x = 3 for down-type and x = 4 for up-type fermions, W the 2x2 sfermion rotation 
matrix, 



hp 



V2 



(e f - If )tasx9 w Z pl + ifZ, 



p2 



f f Rp = -V2e f t&n6 w Z pl . 
If denotes the SU(2)^ isospin and Cf the charge of the fermion / 
The H° k ft ft couplings @ are 

G% x = G^u*^) = (R & Gl R1 (iTf),. = 

,fi / ~7^^u L ~ e uSw) S a+/3 + mw U s C a 



-gR y 



a (^uCq -|- /XSq,) 



+ (l&a) 



(A.8) 
(A.9) 



-rnz e s 2 _ I I ( R ) > (A. 10) 



■)d\T 



G iji = G ( h ° d * dj) = (i? d G d LRl {R 7 
(i| L - e d s 2 w )s a+p + i^Vs 



(/if 



d I cw 



{A d S a + fjLCa) 



^^(A d s a + fic a ) ^e d s 2 w s a+ p + 



(R 



d\T 



2mwCf3 



m w Cp 01 



(A.ll) 



GL = G(H° f* fj) = G(h°f:fj) with a -> a - f 



2 • 



(i. e. sin a = s a — * — c a , cos a; = c a — ► s a , and sin(a + /3) = Sa, +i g — > — c a+ p) 



(A. 12) 



m u (A u cot/3 + //) 



'J 



(A.13) 



G 



1 



m (i (A (i tan/3 + /x) 



(A.14) 



G&4 = G(G° ft fj) = G(A°f*f 3 ) with - - f 



(A.15) 
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(i. e. tan/5 <-> — cot/5). The superscript "u" ("d") denotes an up-type (down-type) 
sfermion. 

The H® Hi f* fj interaction is given by 



with 



rO ttO t-* 



n i 



(A.16) 



ki 



( cos 2 a 
cos a sin a 


V 



cos a sin a 



sin 2 a 











cos 2 (3 
cos P sin (3 








cos (3 sin /? 
sin 2 /3 J 



(A.17) 



A7 



A.:/ 



/ sin a 
—cos a sin a 


V o 



—cos a sin a 
cos 2 a 








sin 2 (3 
-cos (3 sin /3 








\ 



-cos (3 sin (3 

cos 2 /? y 



(A.18) 



dkl 

el 



C kl C kl > 



2 C 2 ^ 



{If - e f s 2 w ) Rl R* + e f s 2 w R{ 2 R 



J2 



(A.19) 
(A.20) 



B Counter term Se 



When we give the renormalized electric charge in the Thomson limit with the measured 
fine structure constant a = el/^n) = 1/137, the counter term 5eo is given by the general 



form [111: 



Sep 
e 



n^(o) 



Sw n^ z (o) 



c w 



(b.i; 



with the momentum derivative of the transverse photon self-energy Il^ A and the 7 — Z° 
mixing self-energy, both for the on-shell photon (p 2 = 0). Fermions and sfermions do not 
contribute to U^ z (0) as a consequence of the fact that the physical photon is massless to 
all orders. However, the contribution of light hadrons to 11^(0) has a large theoretical 
uncertainty p4| |16|j . To avoid this problem, in this work we use the MS running coupling 

(m^)/(47r) = 1/129 as input. The counter term 5e(mz) then 

Q 2 



at Q = m z , <y(m z ) 
becomes 

5e 
e 



2 „2 



e e 



with xj 



(4tt) 2 

rriz for all my < rriz and x t 



A + log ■ 



ml 



(B.2) 



m t . Here A denotes the UV divergence factor. 
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C Proof of the £ independence 



We investigate the wave-function corrections to the process 

A (p)-,x m (h) + x°n(h). (C.i) 

Both the contributions of the transitions A — > G° — > Xm + Xn and A — > Z° — > Xm + X° n 
have a dependence on the gauge parameter £ = £z in the propagators of We 
show that the sum of these contributions is independent of £. 
We start from the matrix elements in a general R% gauge, 

m g = (m AG ) - 1 z (-9Kn4Hhh 5 v(h) , (c.2) 

V ~~ K m z 



v ~ ni z \ P ~ \, ul z J c w 



(C.3) 



The self-energies Hag an d ^-az by (s)fermion one-loop contributions are £ independent. 
The Z x m X°n couplings 0"^ n are 

O-rnn = ~^\~ % raiZ nZ + Z m4 Z n i) . (C.4) 

As limiting cases, A4 G = in the physical unitary gauge £ — > oo and .M z = in the 
£ = (Landau) gauge. Note that the tadpole contributions have to be included |18|, [19 
in II AG . 



We can write MP directly as 

1 

p 2 — £m| 

For M. z we first contract the Lorentz indices 



M G = - 2 gn AG ( P 2 )F^ ni u(hhMk2). (C.5) 



p" (-9,u + (1 - t^T t ) HhWlMh) = f (1 2 - l) u(hWv(k 2 ) , 

V P ~ K m z ) \P ~ K m z ) 

and use u{k\)fy 3 v{k<i) = u{ki){Jfz l + ^ 2 )7 5 t>(&; 2 ) = (m m + m n )u{ki)^ b v{k2)- So we get 



M z = T \ s /^ 2 - 1 — TIaz(p 2 ) 0± (m m + m n ) u(k 1 )<fv(k 2 ) . (C.6) 

- m z \p 2 - £m 2 z J c w 



We use the Slavnov- Taylor identity (see also eq. (3.7)) 

p 2 U AZ (p 2 ) +tm z U AG (p 2 ) = 0, (C.7) 
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and split the sum M° ' + M. z in an obviously £ independent and possibly dependent part, 

M G+Z = ( % 2 ^-U AZ (p 2 )0'^ n (m m + m n ) 
L — mf cw 

+ igU Az ip 2 ) _ p2 . 2 - ^ ~ C) 2 (m m + } ^(feihMfe) ■ (C.8) 

— £m z \ mz p — m z cw J i 

With the relation (proved later) 

m z c w = -0'm n (m m + m n ) , (C.9) 

we get for the part written in the brackets in eq. (|C.8|) 

- (m m + m n )^ (\ + V) =-(rn m + m n ) °/ n V \ ~ im } , (CIO) 
Cw \1Tlz V — m z ) m z c w P — m Z 

and therefore the final result 

M G+Z = -i n AZ (p 2 ) Ol L n {m m + mn) uih^vih) . (C.ll) 

m| c w 

The £ dependence is completely cancelled in (|(J.11| ). 

Finally, we prove ( |C.9| ). With the abbreviation = (Z mi Z n j + Z m jZ ni ) and knowing 
the entries of the neutralino tree- level mass matrix Y (see e. g. eq. (35) in ||11||), one can 
write F^ n4 as 

2 m Z C W F°„ 4 = ^13^31 + *23^32 ~ Y l4 A 41 - F 2 4^42 • (C.12) 

Next we add and subtract the terms I33/I33 + K13A34 and I34A43 + K44A44. Exploiting the 
fact that Y33 = Y44 = and Y34 = Y43, we get 

2 m z c w F° n4 = J2 ( Y k3A 3k ~ Y k4 A ik ) . (C.13) 

k 

Writing the entries of Y in terms of neutralino masses, Y k j = Y^i m iZi k Zij, and using 
Efc Z ik Z jk = Sij we get 

2 mz cw F^ n4 = m\Zi k \Z\ 3 (Z m3 Z nk + Z mk Z n3 ) — Z^ (Z m4: Z nk + Z mk Z n4 )\ 

k,l 

= mi (5 n iZi 3 Z m3 + S m iZi 3 Z n3 — b n iZi 4 Z m4 — b m iZi^Z n ^) 
1 

= m n Z n3 Z m3 + m m Z m3 Z n3 — m n Z n4 Z m4 — m m Z m4 Z n4 

= {Z m3 Z n3 - Z m4 Z ni ) (m m + m n ) = -2 0'^ n (m m + m n ) . (C.14) 

Therefore, eq. (|C.9|) is proven. 

However, from the Slavnov- Taylor identity one can prove in general that the same can- 
cellation of the gauge dependent parts in G° and Z° propagators occurs for any one-loop 
two-body decay of A . 
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